Consequences of a Cosmic Scalar with Kinetic Coupling to Curvature 
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The classical gravitational theory of a scalar field with a gradient coupling to the Ricci tensor is 
examined. This is a scalar- vector-tensor gravitational theory, but in the case that the coupling is 
weak and the scalar evolves like a quintessence field on cosmological time scales, the field equations 
within the solar system are similar to a vector-tensor theory predicting tightly-constrained preferred- 
frame effects. In the early universe, it is shown that strong coupling effects can damp the evolution 
of the scalar field rolling down a potential to help drive an inflationary epoch. In the absence of a 
potential, the strong coupling effects drive a coasting expansion epoch which ultimately terminates 
J"""""- ' in a sudden singularity. 
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(N ; I. INTRODUCTION 

The boundary between the study of new fields and new gravitation has been blurred in the efforts to understand dark 
energy. The observed cosmic acceleration and missing energy problem have lead to speculation of the existence of a 
new field which either acts like a present-day inflaton or moderates a departure from Einsteinian gravity. Consequently, 
scalar-tensor theories of gravity have enjoyed fresh study as candidates to explain dark energy phenomena. In this 
context, we propose to examine a novel kinetic coupling between a scalar field and curvature through the Ricci tensor. 
In particular, we pursue the consequences of gravitational systems described the by Lagrangian 
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wherein gradients of the scalar modulate the strength of gravity. We are motivated in part to find new mechanisms 
to control the evolution of a cosmic scalar field, for inflationary or quintessence scenarios. This theory introduces a 
^ ' new parameter, a, with dimensions of length-squared. Our aim is to study the dynamical behavior and determine 
0^ \ the range for a with respect to cosmological and solar system constraints. 

A gradient coupling is cause for alarm, as the sign of the coupling term in the Lagrangian is indeterminate. If the 
coupling term is sufficiently large and negative to overcome the canonical kinetic energy, then negative energies may 
result. Instabilities may be avoided if a is sufficiently small or if there occurs a separatrix in the equations of motion 
0^ 1 which prevents the coupling from overcoming the canonical energy term. 

The gravitational Lagrangian for scalar-tensor theories, for comparison, is 

°: £s . T= -ffl*) + I ( v«=- m . (2 ) 

To make a correspondence to the Brans-Dicke theory, we set V = and identify the Brans-Dicke scalar as &bd = 
^^■0 2 , with f(R, 4>) = &bdR [H- Another class of scalar-tensor theories dispenses with the scalar field <f> altogether 
and considers a gravitational Lagrangian consisting of a function of the Ricci scalar only, f(R). (See Refs. @, Q for a 
capsule view of this subfield.) Our theory looks like a simple extension of Brans-Dicke to include a vector interaction. 
Next, a broad class of vector- tensor theories are described by the Lagrangian 
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(3) 



(See Ref. Q| for a summary of the properties of scalar- and vector-tensor theories.) The theory we aim to study |T]) is 
obviously a hybrid of the scalar- and vector-tensor theories, if we identify = V M </> and set a 2 = = a 5 = a 6 = 0. 
In fact, our motivation to study the kinetic coupling to curvature derives in part from recent investigations elsewhere 
of the Einstein-aether theory [J] and scalar-vector-tensor theories Q of gravitation. A kinetic coupling to a cosmic 
scalar helps define a preferred frame, and introduces new scalar and vector degrees of freedom to gravity. Our theory is 
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different from the two above-cited cases, which enforce a constraint (V0) 2 = 1, whereas our scalar obeys a generalized 
Klein-Gordon equation. 

In the following, we study the classical theory described by the Lagrangian £ grav (JT() . In section [ill we give the 
equations of motion and discuss solution methods. In section Hill we make a parametrized post- Newtonian analysis to 
show that, for a cosmic scalar field with weak kinetic coupling, the gravitational field equations are similar to that of 
a vector-tensor theory. In section HVl we analyze the scalar field evolution in the case of strong kinetic coupling during 
the early universe. Note that we use metric signature (H ) and curvature convention R^ vaT = d a T^ T — and 

II. EQUATIONS OF MOTION 

The gravitational field equations are obtained by varying the action 

S = j d 4 X y/^g (£grav + ^matter) (4) 

with respect to the metric c/ M „, whereby 
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Here T^ m > IJ,v is the stress-energy tensor for all other matter and radiation. Because the matter and radiation stress- 
energy is conserved independently of the scalar field, and the Bianchi identities still hold, then the scalar field and 
order-a terms form a conserved system. The divergence 

V ^ T (0)^ + T (a)^ = o leadg t0 the scalar field equation f 

motion 

U(j> + V + a [V^iJVV + 2 J R""V M V„0] = 0. (6) 

The field equations ([5]) contain third-derivatives of <fi, whilst the scalar field equation ([6]) contains third-derivatives of 
the metric. 

At this point it is useful to identify a weakly-coupled regime in which 0(a) terms are subdominant to other forms 
of stress-energy, and a strongly-coupled regime in which 0(a) terms dominate. In the weakly-coupled regime, the 
system of equations can be reduced to second-order, by treating the coupling terms perturbatively. When the coupling 
terms and the canonical scalar field stress-energy are subdominant to the matter and radiation, ||T^' W '||, ||TW'*''| <C 
the scalar field equation of motion becomes 

U<j> + V « SirGa [v„S {m) V^ + 2S {m) » u V M V . (7) 

Here we define 

Likewise, we can replace the occurrences of the Ricci tensor in by 8irGS^ m ^ v . As a result, the field equations 

contain at most second-derivatives of the metric. These simplifications are well-suited for the case of a cosmic scalar 
field evolving homogeneously on a cosmic time-scale, as we turn to study gravitation in the solar system. 

III. PPN ANALYSIS 

Given the preponderance of experimental data verifying Newtonian gravity in the weak- field limit, we analyze the 
gravitational field equations to first order beyond the Newtonian approximation, 
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The metric is decomposed as g M „ = r/^ + h^ where rj^v is the Minkowski metric, /i M „ is a perturbation which is at least 
first order in (v/c) 2 , and v is the characteristic velocity of bodies in the solar sytem. Perturbations to the scalar field, 
generated by the stress-energy of the solar system, are written as 8<j>. We will determine the order of 5<j> by comparing 
its field equations to h^. We adopt the standard conventions of the parameterized post-Newtonian (PPN) analysis 
of gravitation (see Ref. [l[), but again note our metric signature. We use the nomenclature 0(1) to indicate <~ (v/c) 2 , 
C(1.5) ~ (v/c) 3 etc. Hence, we expand the field equations (|5I6|) to linear order in /i M „ and <5</>, making assumptions 
as follows. First, (CI) we assume that <f> is spatially-homogeneous and evolves on cosmological time-scales which are 
much longer than the characteristic time-scales of the solar system or laboratory. Consequently, we discard spatial 
derivatives and second and higher time derivatives of <f>. Second, (C2) we assume that the background scalar field is 
a subdominant source of stress-energy in the solar system, whereby HT^^H, | \T i - a ^' J \ < | |T (m)Ail/ | |. Although this 
is inherently a scalar-vector-tensor theory of gravity, we will soon find that these assumptions render it very similar 
to a vector-tensor theory. 

The perturbed field equations are 
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The perturbed scalar field equation of motion becomes 
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And we use the gauge conditions 
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To solve for hoo and hij to 0(1), we collect all of the 0(1) terms in equation (fTU|) . From our constraint (C2) (which 
implies cup 2 < 0(1)) and equation (fTTjl . we see that 8<f> ~ 0(1.5), and thus does not enter into this first order expansion. 
Furthermore, all of the terms proportional to a are also of too high order (being proportional to etc/) 2 ~ 0(1) times 
another perturbative term of order at least one). Therefore, to 0(1), the 0-0 and i-j components of equation (fTU|) are 
unchanged from the case of General Relativity, giving 
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Therefore, the Newtonian limit is preserved. 

The off-diagonal metric components h$i are higher order, 0(1.5), which now includes contributions from the scalar 
field fluctuations, <D(8<t>) Q. Hence, we must solve equation (jlip to leading order before solving the 0-i components 
of equation (fTOf . We rewrite equation (fTTj) as 



□50 = n^d^sr'p + 0(> 2) 
v 2 ^ = <Mr° - i>8T% + 0(> 2) 



V 2 5(f) + 0(2.5) 
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<P (d a h Q0 + dohu) + 4>dih 0i + 0(> 2) 



(14) 



where terms proportional to V and its derivatives have been absorbed into 0(> 2) by assumption (C2), i.e. if is 
subdominant, then V ~ 0(> 1) and 



V"8U 2 + V (-h m tf + 2#0) + V (U8<j> - rTittfJ) ~ 0.{> 2) 



(15) 



Using the superpotential x ( see equation 4.29 of Ref. [l|: V 2 x = 2U) and solutions (fl~3"|) , we find to leading order 

V 2 S<j) = -2^o V 2 X + fahoi. (16) 
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Integrating once gives 

dM = -2j>d 0iX + 4>h 0i ~ 0(1.5) (17) 
Collecting terms of order 1.5 or less in equation (fTO)) . we have for the 0-i component 



^ = SS { Q f + (-2d 0iX + hoi) + 0(2) 

= -pv l + <p 2 (-2doa + h 0t ) (18) 

By our assumption (C2), and the observation that V~ 2 p ~ 0(1), we find that the second term on the right hand 
side of equation (JT5J) will result in a contribution of 0(2.5) to hoi, so that, to post-Newtonian order, hot = |V$ + ^Wi 
where Vi and Wi are defined in equation (4.32) of Ref. [l[. This is the same result derived in General Relativity. 

Finally, we wish to solve for hoo again, only this time, expanding (JTUJ) out to 0(2). Using the stress-energy for a 
generalized fluid, following [l|, then 

^"-/.(i + f + ^ + f + lf 

' 'V 2 U ~ 2V 2 $i - 2V 2 $ 2 - V 2 $ 3 - 3V 2 $ 4 ) (19) 
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where II is the ratio of proper energy to rest energy, p is the pressure, and the potentials $i are defined in equation 
(4.35) of Ref. 0. Next, we expand the Ricci curvature to second order, 



SR 00 (h 2 ) = 5R 00 (h) + ~ 



h l3 dijhoo - dihoodihoo + djh Q0 I dihij + -d h o 



2 

= 6R 00 - V 2 U 2 -4V 2 $ 2 , (20) 
making use of the expressions in the appendix of Ref. [6|] . We find that the zero-zero component of (flTJ|) is 

1 ;V 2 /i 00 - —7} (V 2 C/ 2 + 4V 2 $ 2 ) = -3— (V 2 C/-2V 2 $i-2V 2 $2- V 2 $ 3 -3V 2 $ 4 ) +2a</. 2 V 2 [/ 
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2Z7(1 + 167rGa^) + 2J7 2 - 4$ a + 4$ 2 - 2$ 3 - 6$ 4 (21) 



We absorb the correction to the 2U term into a rescaling of Newton's constant so that, hereafter, the G appearing in 
the PPN potentials (x, Vi, W i: is the G in equation |l|) multiplied by 1 + X = 1 + 167rGa0 2 . Now 



*- = 5 (ifx) + 5 (TtW - + ^ ) (1 ~ X) (23) 

2J7 2 + 4$ 2 1 

(i + x) 2 + TTx 

2U + (2U 2 + 4$ 2 )(1 - 2X) + (-4$i - 2$ 3 - 6$ 4 )(1 - X) (24) 



h °° = 2U+ n ; ' + (-4$x - 2$ 3 - 6$ 4 ) 



where we have again ignored terms proportional to V on the basis of (C2). In terms of the standard PPN parameters, 
we have 

7 = 1 - X P = l-2X £ = 
ai = —6X o?2 = —X Q!3 = — 2X 

Ci = C2 = -5A- Cs = -X (4 = 0, (25) 
as compared to the results for General Relativity 

7-1 0=1 ^=0 
ai = «2 = «3 = 

Ci = C2 = Cs = C4 = 0. (26) 

We note that the scalar field terms involving the canonical kinetic energy and the scalar field potential have been 
assumed or shown to be negligible, contributing at higher orders, so that the theory is essentially General Relativity 
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plus a vector-tensor interaction. However, the PPN parameters derived above (|25|) cannot be easily obtained from the 
generalized vector-tensor theories analyzed in the literature: in most cases a constraint equation from the conservation 
of the field equations (equations 5.50, 5.58 in Ref. Q), has been employed which is inequivalent to the modified Klein- 
Gordon equation (|6llip which we use. The vector-tensor Lagrangian ^ does not contain the kinetic term for the 
scalar field used in our theory {T]). Nevertheless, we obtain effects which are characteristic of vector-tensor theories: 
modification of the degree to which space is curved by mass (7), modification of the non-linearity of gravitational 
superposition (/?), and preferred-frame effects (ai_3). 

A recent survey of American and Russian radar measurements of planetary motion constrains (3— 1 = (0± 1) x 1CP 4 
Q. Similarly, radar ranging of the Cassini probe constrains 7 — 1 = (2.1 ± 2.3) x 10~ 5 10]. These constraints are 
consistent with our assumption (C2). However, they are not the tightest constraint on our theory. By predicting 
non-zero a*, our PPN analysis shows that our gravitational Lagrangian C grav ([1]) specifies a preferred reference frame 
defined by the four- vector V M (/>. Non-zero 03 also threatens non-conservation of momentum [Til ] . By considering 
the orbits of millisecond pulsars in our galaxy, Ref. [ll[ finds the constraint | CK3 1 < 2.2 x 10~ 20 to 90% confidence. 
Consequently, the predicted PPN parameters in this theory can differ from General Relativity by terms which are at 
most of order 10~ 20 , thereby posing a serious challenge to the theory C grav JT]). 



IV. COSMOLOGICAL EVOLUTION 

The equation of motion for the cosmic scalar field <fi in a spatially-flat Robertson- Walker metric is 

• a , ( -a ■ d ■ a a\ 
<j> + 3<f>- + V' - 2a 30- + 3<p— + 6<f> — =0. (27) 



a \ a a a a 

The energy density and pressure are 

P4> = T a + T o = o<P + V + 6a ,, 
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V4> 

However, this system is not universally stable. If we describe the net stress-energy as a perfect fluid with an equation 
of state u>, then small fluctuations obey the equation 

Cl S4> + 3c 2 HS<j> - ^|V 2 (50 + V"6<p = 0, (29) 

where the coefficients are 

ci = l + 3ai? 2 (l + 3w), c 2 = l-3aH 2 w(l + 3w), c 3 = 1 - 3aH 2 {l - w). (30) 

In the weak-coupling regime, with ^\a\ ^ H~ x , the system is stable against the rapid growth of small perturbations. 
In the strong-coupling regime, with y/\a\ 3> H ~ 1 , the system is unstable for w > — 1/3 because c\ and C3 have opposite 
signs. Conversely, stability during the standard radiation- and matter-dominated epochs requires y/\a\ < H^ 1 
throughout. The most constraining bound occurs at the beginning of the radiation epoch. When w < —1/3 as occurs 
during inflation, c\ and c 3 have the same sign which ensures stability, regardless of the magnitude of a. In fact, we 
note that the coefficient c\ multiplying the second derivative of either ip or Sip changes sign across w = —1/3 in the 
strong-coupling limit, suggesting the possible existence of a separatrix, marking a boundary across which the scalar 
field cannot evolve. 

Given the stability criteria, we are motivated to consider the role of the kinetically-coupled scalar field during 
inflation. We consider a scalar field potential V = \m 2 4> 2 in an inflationary background, with scale-factor a oc 
exp (Ht) . In this case, the (/(-equation becomes 



+ 3ff0+ „ 2 = O (31) 

1 — oaH z 

so that the kinetic-coupling renormalizes the mass. Defining a = —aH 2 , then in the limit a 1 (LI) the effective 
mass is dramatically reduced, m — > m/\/6a. Hence, the kinetic drag terms observed in equation (|2~T|) slow the field 
evolution in a way identical to a flattening of the potential. Taking the limit a 3> m 2 / H 2 (L2) then the two solutions 
to the scalar field evolution are <p oc exp (Xt) with Ai = — 3H and A2 = —m 2 /18aH. Whereas the first solution decays 



6 



rapidly, the second solution describes the slow-roll of the scalar field down the potential. As well, the limits (Ll-2) 
ensure that the energy density and pressure are dominated by the scalar field potential, so that ~ — p^. Hence, a 
self-consistent inflationary solution is possible. 

In the absence of a potential, setting V = 0, the system also has interesting solutions. We further neglect any 
background radiation or matter, and recast the field equations as 

H' _ v' 1 1 
IF ~~ ~ ~ 2av 2 + \2aH 2 

"" = -I (»' 2 + ' 4 + § > w + < 3 + 4>" 2 - jsp - i i3 + 4>) < 32 > 

where v = vo7rG</>, x = In a, and ' = d/dx. We focus on the strong-coupling regime with a < and start the field 
with \a\ 3> H 2 3> v 2 , v > 0, and w < —1/3. Following the scalar field dynamics numerically, we find that the field 
approaches the analytic solution 



ur M/^ + ^m a<V2oC, (33) 
V 2a a A 

where C is determined by the initial conditions. This analytic solution corresponds to a coasting universe with 
w = —1/3 and a oc t. The full system, however, never reaches this ideal solution. As a — > a c = \/2aC, the full 
solution grows singular: v vanishes like oc y^log a c /a whereas v' and - diverge like oc 1/ \/log a c /a. This indicates a 
sudden singularity [T2l . [l3| of the curvature. As the critical value of the scale factor is approached, in finite time the 
energy density remains finite and positive but the negative pressure grows divergent, corresponding to an equation of 
state w — > — oo. It seems clear that the kinetic coupling leads to a dangerous instability. 

V. DISCUSSION 

The consequences of a kinetic coupling between a cosmic scalar field and the Ricci curvature tensor have been 
studied in this work. In the weak coupling regime, the theory is severely constrained by precision tests of General 
Relativity Solar system measurements limit the contribution by the coupling to the energy density and pressure of a 
cosmic scalar field to be smaller than ~ 10~ 6 of the critical density, which eliminates a possible role in a quintessence 
scenario. Although this theory is a hybrid of scalar- and vector-tensor theories, we find that the theory suffers from 
the same problems as vector-tensor theories. Hence, the coupling stress-energy density is restricted to be smaller than 
~ 10~ 20 due to limits on preferred-frame effects. 

A possible new mechanism to control the rolling of a cosmic scalar field has been our primary motivation to 
investigate this theory. Indeed, we have found that in the strong coupling regime, with large and negative a, during 
an inflationary epoch, the coupling terms act like a brake on the scalar field evolution and effectively flatten the slope 
of the scalar field potential. This behavior may suggest a possible role for such a coupling in the early universe. 

The theory also displays dangerous instabilities. In the absence of a potential, in the strong coupling regime, with 
large and negative a, the coupling leads to a sudden singularity in the curvature. The strong coupling regime is 
non-viable in standard radiation- or matter-dominated epochs due to an instability to fluctuation growth. These 
pathologies suggest a strong kinetic coupling may be unsuitable for a healthy, long-lived universe. 
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